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Abstract 

We study variations of the lifetimes of high-l 7 solar p modes in the quiet and 
active Sun with the solar activity cycle. The lifetimes in the degree range £ = 
300-600 and v = 2.5-4.5 mHz were computed from SOHO/MDI data in an area 
including active regions and quiet Sun using the time - distance technique. We 
applied our analysis to the data in four different phases of solar activity: 1996 (at 
minimum), 1998 (rising phase). 2000 (at maximum), and 2003 (declining phase). 
The results from the area with active regions show that the lifetime decreases as 
activity increases. The maximal lifetime variations are between solar minimum 
in 1996 and maximum in 2000; the relative variation averaged over all £ values 
and frequencies is a decrease of about 13%. The lifetime reductions relative to 
1996 are about 7% in 1998 and about 10% in 2003. The lifetime computed in the 
quiet region still decreases with solar activity although the decrease is smaller. 
On average, relative to 1996, the lifetime decrease is about 4% in 1998, 10% in 
2000, and 8% in 2003. Thus, measured lifetime increases when regions of high 
magnetic activity are avoided. Moreover, the lifetime computed in quiet regions 
also shows variations with the activity cycle. 
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1. Introduction 

The excitation of solar acoustic oscillations is attributed to turbulent convection 
beneath the solar surface. The structural and dynamical properties of the Sun, 
as well as the processes of excitation and damping of the acoustic modes, change 
with the solar cycle. Variation in the excitation and damping are observed as 
changes in the mode amplitudes and lifetimes. The study of variations in these 
parameters provides information on the mechanism driving the solar oscillations. 
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The p-mode lifetime in global helioseismology is computed from the line width 
as (2ttj)^ 1 , where 7 is the FWHM of the p-mode line profile in power spectra. 
Several authors flChaplin et al., 2000[| Jimenez- Reyes et al., 2003[|Howe et al., 2003| )[ 
have shown significant temporal variations in the mode widths using low-£ data. 
They reported a 23% increase in the mode widths (i.e. decrease in the mode life- 
times) with increasing solar activity. They also found that the variations are fre- 
quency dependent, being largest near 3.1 mHz, and independent of £. The results 
have been confirmed with intermediate-^ data, showing a 14% increase in the 
widths with activity (Komm, Howe, and Hill, 2000; Salab ert and Jimenez-Reyes, 2006] ) .[ 

The authors also analyzed changes in the mode amplitudes, power, and energy-| 
supply rate from solar minimum to maximum and concluded that the varia- 
tions probably arise from an increase of the damping only, since the net forc- 
ing of the modes remains constant. Some of the authors noted that, even if 
no significant changes in global p-modc excitation are evident, it may possi- 
bly change locally. Possible mechanisms, in addition to a continuous turbulent 
source of excitation, are transient phenomena such as flares (Haber et al., 1988; 



Kosovichev and Zharkova, 1998 Donea, Braun, and Lindsey, 1999||Ambastha, Basu, and Antia, 2003| | 



Martinez-Oliveros, Moradi, and Donea, 2008l|Karoff and Kjeldsen, 2008[ ) and CMEs| 
( Gavryusev and Gavryuseva, 1999 ). 

One potential source of damping of the acoustic waves is magnetic activity 
that is present on the solar surface. Results of global ( Komm, Howe, and Hill, 2002[ )[ 



and local (Rajaguru, Basu, and Antia, 2001 Howe et al, 2004) analyses have 



shown a strong dependence of the mode amplitudes and lifetimes on the loca- 
tions of active regions, with amplitude and lifetime decreasing in the presence of 
strong magnetic fields. Also, in contrast to the global analysis, Howe et al. (2004) 
reported a significant decrease of the energy supply rate of p modes in active re- 
gions. It is known that the power of the oscillations in sunspots and plages is two- 
three times lower than in the quiet Sun ( |Woods and Cram, 1981| ). Possible mech- 
anisms that might be responsible are the suppression of excitation in sunspots as 



a strong magnetic field inhibits convection (Parchevsky and Kosovichev, 2007), 



absorption of the p-mode waves in sunspots (Braun, Duvall, and LaBonte, 1987 



|Chen et al., 1996} |Chou et al., 1996| , the different height of spectral-line forma- 
tion due to the Wilson depression, and the modification of p-mode eigenfunctions 



by the magnetic field (Hindman, Jain, and Zweibel, 1997). 

Another possible contribution to p-mode damping is variations in the con- 
vective properties near the solar surface. Most likely, these variations arise from 
the influence of magnetic structures, as thermal perturbations in the convection 
zone cannot alone explain observed changes in irradiance and p-modc frequencies 
over the activity cycle ( jBalmforth, Gough, and M crryfield, 1996). Houdek et al. 
(2001) made an attempt to connect changes in the horizontal length scale of 
convective eddies and the p-mode damping rates and found that decreasing the 
horizontal size of convective eddies increases the p-mode damping. 

Whatever mechanisms generate and damp solar p modes and drive the solar 
cycle, they are believed to operate in the convection zone. High-degree (£ > 200) 
acoustic waves propagate in the upper convection zone and yield information 
about inhomogeneitics in structure and dynamics of this region. Determination 
of the line widths and other parameters of the high-£ p modes has been a problem 
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because of the blending of different modes in power spectra. Burtseva et al. 
(2007) compared the lifetime of high-^ p modes at solar minimum and maximum 
using an alternative method based on time - distance analysis first proposed by 
Chou et al. (2001). The conclusion of Burtseva et al. (2007) is consistent with 
low- and intermediate-^ analysis (i.e., the lifetime decreases with increasing solar 
activity). In this work we use the same method to measure lifetimes of high- 
degree solar p modes in the range £ = 300 — 600 and v = 2.5 — 4.5 mHz in four 
different phases of solar cycle in an area including active regions and a quiet-Sun 
area and discuss the results. 



2. Method of Lifetime Measurement 

The lifetime measured with time - distance analysis is the lifetime of a wave 
packet, not of an individual mode. The wave-packet lifetime is determined by 
changes in the amplitude and width of the cross-correlation function of the packet 
with time. The amplitude of the cross-correlation function decreases exponen- 
tially with the number of skips of the wave below the photosphere. This effect has 
been interpreted as the dissipation of solar p-mode power flChou et al., 200lj ). 
As a wave packet consists of many p modes over some range of phase velocities 
(2iti>/£) dispersion causes it to extend (i.e., its amplitude decreases and its width 
increases). In previous studies ( |Chou and Ladenkov, 2007"||Burtseva et al., 2007[ ),| 
it was found that the width of the cross-correlation function increases with 
the number of skips, demonstrating the effect of dispersion of the wave packet. 
Dispersion does not change the energy of the wave packet, which is the product 



of amplitude squared and width of the wave packet (Jackson, 1975). Taking 
into account both the dissipation and the dispersion, we can define the life- 
time as the e-folding time of the energy of the wave packet due to dissipation 
QChou and Ladenkov, 2007||Burtseva et al., 2007[ ): 

A 2 n W n = A 2 W e- n ^/ T , (1) 

where A n and W n are the amplitude and the width of the cross-correlation of 
skip n, respectively, r en is the one-skip envelope time, and T is the lifetime of 
the wave packet. 



3. Data Analysis 

In this study we analyzed two 512-minute scries of MDI Dopplcrgrams in each 
of the following epoch: at solar minimum (24 May and 18 September 1996), in 
the rising phase of activity (5 and 6 February 1998), at solar maximum (9 and 
11 September 2000), and in the declining phase of activity (12 and 13 December 
2003). 

We obtained information on the distribution of magnetic features and the 
magnetic-field strength on the analyzed days using MDI 96-minutc magnetograms| 
and acoustic-power maps. The areas including active regions and the quiet areas 
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Table 1. Unsigned magnetic field (in Gauss) averaged over the whole area 
including active regions and the quiet area (see Figure 1). 



Date 


Area including 


Quiet Area 


Area including 


Quiet Area 




active regions 


(MDI) 


active regions 


(KPM) 




(MDI) 




(KPM) 




24 May 1996 






1.73 ± 0.03 


2.49 ± 0.11 


18 Sep 1996 


2.43 ± 0.02 


1.90 ± 0.07 


1.98 ± 0.03 


2.80 ± 0.11 


05 Feb 1998 


3.60 ± 0.04 


1.71 ± 0.04 


3.25 ± 0.07 


2.15 ± 0.11 


06 Feb 1998 


4.45 ± 0.06 


2.02 ± 0.07 


3.01 ± 0.08 


2.05 ± 0.13 


09 Sep 2000 


15.03 ± 0.33 


5.40 ± 0.12 


15.00 ± 0.12 


6.75 ± 0.20 


11 Sep 2000 


16.21 ± 0.33 


4.93 ± 0.17 


15.20 ± 0.11 


5.42 ± 0.16 


12 Dec 2003 


3.99 ± 0.26 


3.37 ± 0.14 






13 Dec 2003 


5.85 ± 0.27 


3.53 ± 0.14 







analyzed in this work are shown in Figure 1. Concentrations of active regions 
are denoted with black contours, which correspond to the magnetic field level 
of 25 G after smoothing of the images. It is seen that there is no concentration 
of active regions in the quiet areas. On both days in 1996 there were no active 
regions present on the solar disk. On the days in 2000 the Sun was active with 
extended active regions, with a relatively low maximal value of magnetic field of 
800 G. On the days in 1998 and 2003 the central part of the Sun was quiet, but 
a few small sunspots with magnetic field up to 2000 G were present on the disk 
at about 40° from the disk center. 

We characterize the activity level on the analyzed days by computing a mag- 
netic activity index (MAI). The lcr noise level in a single MDI magnetogram 



is about 20 G (Scherrer et al, 1995). To reliably compute MAI we first set all 



unsigned flux less than 2.5<r = 50 G to zero and correct images for outliers 



using a technique similar to that described in (Basu, Antia, and Bogart, 2004). 
We compute MAI by summing the absolute values of the magnetic flux inside 
the analyzed areas and normalizing the sum by the total number of pixels in 
the areas instead of the number of nonzero pixels. We also compute MAI using 
NSO/Kitt Peak magnetograms (KPM) for comparison. The lcr noise level of 



KPM data is about 5 G ( Wenzler et al, 2004[ ). Similarly, we set all unsigned 



flux less than 2.5ct = 12 G to zero. Then we multiply the KPM images by 1.7 
to correct for the difference in scale as MDI shows a larger magnetic flux than 
KPM. The computed MAIs are shown in Table 1. The MDI magnetograms for 
24 May 1996 and KPM data for the two days in 2003 are absent. 

The data analysis and cross-correlation function computation procedure is 
similar to that used in ( |Burtseva et al, 2007[ ) and briefly described as follows: 

i) Images are remapped, tracked, filtered with 15-minute running mean, and 
transformed into the (£, m, v) domain, where m is the azimuthal degree. 

ii) A Gaussian filter of FWHM = 2.0 mHz centered at a frequency v$ is applied. 
in) A phase-velocity filter is applied to isolate the modes in a range of angular 

phase velocity (w = 2ttv/£) forming a wave packet with the central fre- 
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Figure 1. MDI magnetograms of the analyzed days. Lifetimes were computed using the whole 
area (99.4° X 85.2°), including active regions, and in the quiet area (45.9° X 15.7°) shown as 
a white box around the disk center. Black contours denote active-region concentrations. 

burtseva09_4.tex; 21/06/2009; 1:40; p. 5 



6 



O. Burtseva et al. 



quency v$ and the corresponding degree £$. The filter is smoothed with a 
Hanning window in frequency. 

iv) The data are reconstructed back to the space-time domain. 

v) The center-annuli cross-correlation C(A,r) is computed with 

C(A,r) = y*(0,t)*(A,t + T)dt, (2) 

where ^(O,^ is the signal measured at the central point at time t and 
^(A, t + r) is the signal averaged over an annulus at an angular distance A 
from the central point at time t + T. The procedure is repeated for different 
central points at the disk center. The cross-correlation functions of different 
central points are averaged. 

vi) For each skip, C(A, r) is fitted by a Gabor wavelet ( Kosovichev and Duvall, 1997[ )[ 



G(t) — A ■ cos[27w(t — T p h)\ ■ exp 



(T - T, 



\2 



en } 



2W 2 



(3) 



where A, v, and W are the amplitude, frequency, and width of a Gaussian 
envelope, respectively. The parameters r p h and T en are the phase time and 
envelope time, respectively. For each skip, A, W, and r cn are averaged over 
five angular distances near the maximum of the cross-correlation amplitude 
to reduce noise. The procedure is repeated for both positive and negative r. 
The values of A and W are the averages of results for positive and negative 
r, since they are found to be very close in magnitude. 
vii) The lifetime (T) is determined by fitting In (A 2 1F) versus n according to 
Equation (1). The details are discussed in Burtseva et al. (2007). 

The above analysis was performed for the whole area including active regions 
and the quiet area near the disk center shown in Figure 1. A minimum four- 
skip distance is required to compute lifetimes with this method. The size of the 
quiet area does not allow us to perform the analysis for the wave packets with 
Iq = 300, uq = 2.5 — 4.5 mHz and £q = 400, v$ = 4.0 and 4.5 mHz as the distance 
required for the lifetime computation for those waves is from about 13° to 24°. 
For the quiet area we have about 7° maximum, thus we are able to compute 
lifetimes only for the wave packets £q = 400, uq = 2.5 — 3.5 mHz £q = 500 — 600, 

= 2.5 — 4.5 mHz which requires distances from about 5° to 7°. 



4. Results and Discussion 

The lifetimes of solar p modes computed in the area including active regions are 
shown in Figure 2. It is seen that the lifetime decreases as activity increases. The 
results are in agreement with all previous studies of the p-modc lifetimes men- 
tioned earlier. The percentage of the lifetime decrease for all analyzed days 
relative to the days in 1996 is listed in Table 2. The maximal lifetime variations 
are between solar minimum in 1996 and maximum in 2000; the relative variation 
averaged over all £ values and frequencies is about 13%. The lifetimes in 2003 
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Figure 2. p-mode lifetime computed in the area including active regions, averaged over the 
results from two 512-minute time series in each epoch versus frequency for various Iq. The 
frequency is determined from the Gabor fit for the first-skip cross-correlation function. 
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are closer to the lifetimes in 2000, and the lifetimes in 1998 arc closer to those 
in 1996. The lifetime variations relative to 1996 are about 7% in 1998 and about 
10% in 2003. As seen in Table 1, the activity level was higher on the analyzed 
days in 2003 than in 1998. 

The lifetimes computed in the quiet area are presented in Figure 3. After 
excluding active regions from the analysis, the lifetimes decrease with solar 
activity, although the decrease is smaller. On average, relative to 1996, the 
decrease in the lifetime is about 4% in 1998, 10% in 2000, and 8% in 2003. 
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The daily variations of the lifetime in 1996 are within the uncertainties of the 
measurements. 

Thus, excluding high magnetic activity regions increases the measured life- 
time. From the difference between lifetimes and the magnetic field strength of the 
quiet area and the area including active regions, and also from the distribution 
of the active features on the solar disk on the analyzed days (see Figure 1), we 
speculate that acoustic waves are affected more by large active regions as seen 
at solar maximum, and less by small sunspots with strong magnetic field. This 
is likely because the wavelengths of the acoustic waves are typically larger than 
the smaller sunspots and better matched with the large active regions. 

The lifetime computed in the quiet area still shows variations with activity 
cycle. What could cause the variations in the quiet area where there are no large- 
scale magnetic field concentrations? One explanation could be that the convec- 
tive properties change with solar cycle. As discussed in the Introduction, Houdck 
et al. (2001) found that damping of p modes increases as horizontal size of 



convectivc eddies decreases. It has been reported by some authors (Muller, 1988 



Berrilli et al., 1999) that the horizontal scale of the solar granules decreases from 



solar minimum to maximum, which perhaps indicates that the horizontal size 
of the convective eddies docs indeed decrease. Meunier et al. (2008) reported 
that supcrgranular cells are smaller at solar maximum. However, other studies 
arrive at the contradictory conclusion that granular size increases with activity 
(Hanslmeier and Muller, 2002 Saldana-Munoz, Muller, and Hanslmeier, 2004 1. 

Solar-cycle variation of the lifetime in the quiet area caused by purely struc- 
tural perturbations in the convection zone should result in lifetime variations 
with no correlation with magnetic field strength. To test this we computed the 
average lifetime for each wave packet over the four epochs and normalized each 
of the lifetimes by the average. Then we plot the normalized lifetimes of all wave 
packets analyzed versus MAI computed in the area including active regions and 
the quiet area in Figure 4. The lifetimes obviously anticorrelate with MAI for the 
area including active regions. In the quiet area we also see an inverse correlation 
of the lifetimes with MAI. This suggests that the magnetic field somehow plays 
a role in the lifetime variation of p modes in the quiet Sun. Future work will 
involve more statistics to confirm the conclusions. 

All of the analyzed quiet regions have different levels of magnetic field (see 
Table 1). Some isolated pixels in the quiet regions have magnetic fields as strong 
as «500 G, which are the strong-field component of the network field at the 
supergranular boundaries ( |Harvey, 1971|[Lm, 1995[ ). 



Table 2. Percentage of the lifetime decrease relative to 1996, averaged 
over all I values and frequencies of 3.0, 3.5, and 4.0 mHz. 



Year/ 1996 


Area including active regions (%) 


Quiet Area (%) 


1998/1996 


6.9 


4.2 


2000/1996 


13.2 


9.6 


2003/1996 


9.6 


7.6 
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Figure 4. Normalized lifetimes versus magnetic activity index in the area including active 
regions (left) and the quiet area (right). Top and bottom panels correspond to the MAI 
computed from MDI and Kitt Peak magnetograms, respectively. 



Thus, the elements of heightened magnetic activity in the quiet regions could 
cause the residual lifetime variations in the quiet area. The MAIs of the analyzed 
quiet areas also show solar-cycle variations. One can clearly see the increase 
in MAI from Table 1 (plotted in Figure 5) corresponding to the days at solar 
maximum for both the area with active regions and the quiet area. Most authors 
who analyzed solar-cycle variations of the irradiance and photosphcric magnetic 
field reported little overall solar-cycle variations in the quiet-Sun network ele- 
ments arising from decaying active regions, and no variations in the intranetwork 
field ([HarTey, 1994| |Topka, Tarbell, and Title 19971 |Pevtsov and Acton,"200Tj 



|Hagenaar et al., 2003||Sanchez Almeida, 2003||Pauluhn and Solanki, 2003[ ). Muller| 

et al. (2007) reported solar-cycle related variations of the granulation contrast 
at the solar surface connected to the temperature differences between granules 
and intergranular lanes. The contrast variations may be due to changes in the 
physical structure under the surface or a varying amount of magnetic flux in the 
intergranular lanes. 

The recently discovered horizontal magnetic field in quiet regions in the solar 
photosphere ( |Harvey et al., 2007| |Lites et al, 2008D might also affect the prop- 
agation of the acoustic waves in the quiet Sun in a way different from that 
of the vertical component of the magnetic field. As the authors speculate, the 
possible driving mechanism for the horizontal field is granular and supergranular 
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Figure 5. Temporal evolution of MAIs from Table 1 vs. four epochs used in the lifetime study. 
The line and symbol styles arc shown in the top-left corner of the figure. The plot clearly shows 
solar cycle variations of the MAI in both the area including active regions and the quiet area. 



convection and transformation of the field lines in response to the evolving flux 
distribution in the network and intranctwork field in the quiet Sun. 

In future work, we will analyze a large number of quiet areas to confirm 
temporal variations of the lifetimes of high-£ acoustic waves found in this work. 
We will also study variability of the apparent power of the acoustic modes in 
an area including active regions and a quiet area in different phases of the solar 
cycle. 
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